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Abstract 

We discuss properties of conformal geodesies on general, vacuum, and 
warped product space-times and derive a system of conformal devia- 
tion equations. The results are used to show how to construct on the 
Schwarzschild-Kruskal space-time global conformal Gauss coordinates which 
extends smoothly and without degeneracy to future and past null infinity. 
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1 Introduction 



Geometrically defined coordinate systems are convenient in the local analysis of 
fields but they often degenerate on larger domains. It is well known, for instance, 
that Gauss coordinates are of limited use because the underlying geodesies tend 
to develop caustics. In |^] conformal Gauss coordinates have been introduced 
and successfully employed in local studies of conformally rescaled space-times. 
Here the time-like coordinate lines are generated by conformal geodesies, curves 
which are associated with conformal structures in a similar way as geodesies 
are associated with metrics. In this case the danger of a coordinate break-down 
is even more severe, because the conformal geodesies generating the system 
can not only develop envelopes or intersections, as occur in caustics of metric 
geodesies, but they may even become tangent to each other. However, it will 
be shown in this article that the structure responsible for this difficulty also has 
useful aspects. 

Subsequent analyses of conformal Gauss coordinates revealed that conformal 
geodesies possess various remarkable properties. In they have been studied 
systematically in the context of the conformal field equations. Though this at 
first introduces additional complications into the equations because it requires 
the use of Weyl connections, it leads to unexpected simplifications. While the 
Bianchi equation satisfied by the rescaled conformal Weyl tensor, which is given 
by (f jki — C" jki where Q denotes the conformal factor, always implies 
hyperbolic reduced systems of partial differential equations, the generalized 
conformal field equations (in vacuum, possibly with a cosmological constant) 
admit in the new gauge hyperbolic reductions in which the frame coefficients, 
the connection coefficients and the components of the conformal Ricci tensor 
are governed by equations of the form 

9r e'" fc = . . . , drfi^k^---, drRjk^---, 

where r denotes the time variable in the given gauge and the right hand sides 
are given by algebraic functions of the unknowns j, , ^ , Rjk , ju and 
known functions on the solution manifold. 

Moreover, associated with a conformal geodesic is a function along that 
curve, the conformal factor, which is determined up to a constant factor that 
can be fixed on a given initial hypersurface. This function necessarily has zeros 
at points where the conformal geodesic crosses null infinity. It turns out that for 
given data on the initial hypersurface this function can be determined explicitly 
and acquires a form which allows us to prescribe to null infinity by a suitable 
choice of initial data a finite coordinate location. 

In [Q these facts are used to provide a rather complete discussion of Anti-de 
Sitter- type solutions, including a characterization of these solutions in terms of 
initial and boundary data. In these properties form a basic ingredient of 
a detailed investigation of the behaviour of asymptotically fiat solutions near 
space-like and null infinity under certain assumptions on the Caucliy data for 
the solutions. In particular, the cylinder at space-like infinity introduced in 

to remove the conformal singularity at space-like infinity is obtained as a 
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limit set of conformal geodesies. The results obtained there strongly suggest 
that even under the most stringent smoothness assumptions on the conformal 
data near space-like infinity the solutions will in general develop logarithmic 
singularities at null infinity. However, they also suggest that these singularities 
can be avoided if the data satisfy certain regularity conditions at space-like 
infinity. 

The conversion of these conjectures into facts requires the proof of a certain 
existence result for the regular finite initial value problem at space-like infinity 
formulated in While the basic difficulty of this proof has nothing to do 
with the specific features of the underlying conformal Gauss system, only the 
proof will show that the latter is doing what we expect it to do: to extend 
near space-like infinity smoothly and without degeneracy to null infinity if null 
infinity admits a smooth diffcrcntiable structure at all. 

An independent proof that conformal Gauss systems behave this way on gen- 
eral space-times is difficult, because it requires the information on the asymp- 
totic structure of the solution which we first hope to obtain by completing the 
analysis of [^. The question whether conformal Gauss systems can be used 
globally is even more difficult. 

In the context of the hyperboloidal initial value problem with smooth initial 
data it can be shown that these coordinates remain good for a while and that 
they can in fact be globally defined for data sufficiently close to Minkowskian 
hyperboloidal data. However, nothing is known for data which deviate strongly 
from Minkowskian once. There remains the possibility to study the behaviour 
of conformal Gauss systems on specific solutions. It has been shown in ||^ that 
there exist conformal Gauss systems on the Schwarzschild space-time which 
smoothly cover neighbourhoods of space-like infinity including parts of null in- 
finity. 

One may be tempted to think that this is only possible due to the weakness of 
the field in that domain (although it doesn't look weak in the conformal picture) 
while in regimes of strong curvature a break-down of the coordinates will occur. 
It turns out that this is not true. It is the main results of this article that 
there exist conformal Gauss coordinates on the Schwarzschild-Kruskal space- 
time which smoothly cover the whole space-time and which extend smoothly 
(in fact analytically) and without degeneracy through null infinity. 

This raises hopes that conformal Gauss systems stay regular under much 
wider circumstances than expected so far and that we may be able to exploit the 
simplicity of the reduced equations and the fact that conformal Gauss systems 
provide by themselves the conformal compactification in time directions under 
quite general assumptions. How robust these systems really are, whether they 
will stay regular under non-linear perturbations of the Schwarschild-Kruskal 
space, remains to be seen. Though a general analytical investigation of these 
equations appears difficult at present, important information about them may 
be obtained by numerical calculations. Conversely, there are good reasons to 
believe that it should be possible to give strong analytical support to numer- 
ical work based on the use of conformal Gauss systems because the latter are 
amenable to a geometric analysis. 
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Wc begin the article by discussing properties of conformal geodesies on gen- 
eral space-times and then specialize to vacuum space-times, using also some of 
the results which have been obtained already in the articles referred to above. 
For later application wc work out the conformal geodesic equations on warped 
product space-times. Various features of conformal Gauss systems are then 
illustrated by explicit examples on Minkowski space-time. 

After specializing further to warped product vacuum space-times a confor- 
mal Gauss system on the Schwarzschild-Kruskal space-time will be analysed. 
Expressions for the conformal geodesies determined by suitably chosen initial 
data arc derived in terms of elliptic and theta functions. It is shown that the 
curves extend smoothly through null infinity and that the associated conformal 
factor defines a smooth conformal extension if the congruence is regular. We 
do not attempt to derive the explicit expression of the rcscalcd metric in terms 
of the new coordinates, because this cannot be done anyway in more general 
situations. Instead, the regularity of the coordinate system and the conformal 
extension is established by analysing a system of equations which is the ana- 
logue for conformal geodesies of the Jacobi equation for metric geodesies. While 
specific features of the Schwarzschild solution are used here, it is clear that sim- 
ilar techniques apply in more general cases. Wc end the article by indicating 
as a possible application the numerical calculation of entire asymptotically flat 
vacuum solutions. 



2 Conformal geodesies on pseudo-Riemannian 
manifolds 

Before we introduce conformal geodesies we recall a few concepts and formulae 

of conformal geometry. On a pseudo-Ricmannian manifold (M, g) of dimension 
n > 3 we consider two operations preserving the conformal structure defined by 
g: (i) conformal rescalings of the metric 

9tj,v ^ gfiv = 9p.u, (1) 

with smooth conformal factor 51 > 0, (ii) transitions V ^ V of the Levi-Civita 
connection V of a metric g in the conformal class of g into (torsion free) Weyl 
connections V with respect to ^. In terms of the Christoffel symbols F and the 
connection coefficients F, defined by V^^ =Y^'' vdp and Vg^ = ^n^ udp 
respectively, the transition is described by 

F^ ^1/ > F^ 1/ = F|j ^ V -\- S[f)^ P V, (2) 



S[j)^ % = S%U + SP^f^+g^^ gp^ fx, 

with a smooth 1-form /. For given / the difference tensor S{f) depends only on 
the conformal structure of g. The transformation of the Levi-Civita connection 
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f of g under |^ is a special case of in which the 1-form is exact 

r^.%^r^p,^^^% + s{f)/\ with f^n-Un. (3) 

Conversely, if the 1-form / is closed (||) arises locally from a rescaling of g with 
a suitable conformal factor. By (||) we have 

^ p9iiu = -"^ Ip9t,v, (4) 

It follows from this that V preserves the conformal structure of g in the sense 
that with a function 9 > Q satisfying on a given curve x{t) in M 

ViO^e <f,x>, (5) 

the metric O^g^u is parallely transported along x{t) with respect to V. 

The curvature tensors of the connections V and V, defined by (Va Vp — 
Vp Vx)X'' = Rt" ^xp etc., are related by 

R'' .Ap - i?^ .xp = 2 {V[A Sp]^, + Ss ^ [A Sp] ^ J (6) 
were indices are raised or lowered with respect to g. The tensor 
1 . A n - 2 - 1 



Lnu — 



{ R(p^) — R[p,y] - 2(n-l) ^ ^' 



which occurs in the decomposition = 2 {g^^ [>, Lp^y—g'^ ^ L\^xp\— gu[\ Lp] ^}+ 

i/Ap of the curvature tensor of V into its trace parts and the trace-free, con- 
formally invariant conformal Weyl tensor ^\p, is related to the tensor 

by the equation 

Vi. - fpfv+ gp.u -^fxf^^ Lp^ - Lpu. (7) 
2.1 Conformal geodesies 

A conformal geodesic associated with the conformal structure defined by g is 
given by a curve x{t) in and a 1-form &(t) along x(t) such that the equations 

{\JixY + S{b)x\x^±P (8) 



(Vi.6),--5^5(6)A^,i^=iA.i\ (9) 

hold on x{t). For given initial data a;* G M, i* € T^^M, G T*^M there 

exists a unique conformal geodesic x{t)^ ^(t) near x■^, satisfying for given G R 

a;(T*)=x*, i(T*) = i*, b{T^,) = b^,. (10) 
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From (||) it follows that the sign of g{x, i) is preserved along x(t), since we have 

Vi(.g(i, i)) = -2<b,i> g{x, x). (11) 

In particular, if g{x, x) = holds at one point it holds everywhere on x{t) and 
(H) implies that x{t) can be reparametrized to coincide with a null geodesic of 

9- 

Let x{t), &(t) and a;(f), b{f) be two solutions to (||), with g(x,x) ^ 0. 
We derive now the conditions under which the curves x(t)^ x{f) coincide locally 
as point sets such that there exists a local reparameterization t = t(t) with 
x{T{f)) = x{f). The latter relation and (Q) imply xd^r + 2 < b — b,x > 
X (9tt)^ — g{x, x)(b — b) {dfT}'^ = which is equivalent to the equations 

b-b^ax\ d^T + ag{x,x){drTf ^0, (12) 

with some function a. In the first of these equations the index of x is lowered 
with the metric g. From the first equation and follows 

a = 2a < b,x > + — a^ g{x, x), (13) 

which together with (|l^) is equivalent to our requirement. Equations (|ll|), ( p^ 
imply dr{ag{xjx)) = 1/2 {ag{x,x))'^ which has the solutions 

.^ 2a* g{x*,x*) 

ag{x,x) = 



2 — a* g(i*,i*) At' 

where A r = r — and A f = f — . From these solutions follows finally with 
& 

4eAf - 1 2a.j, q(x.j,,x*) 

At= — . , , _ , b = b+——- *\ *' X, 14 

1 + 2 ea, g[x*,x*) Ar g(x,x) 2 — a* gyx^.x^,) Ar 



e, a*, r*, e R, 0. 

Thus, the changes of the initial data (jl^) which locally preserve the point set 
spread out by the curve x{t) are given by 

i*— ^4ei*, 6* ^ 6* + a* i, e,a*GM, e 7^ 0. (15) 

The 1-form remains unchanged and an affine parameter transformation results 
if a* = 0. With suitable choices of the free constants, however, any fractional 
linear transformations of the parameter can be obtained and it can be arranged 
that r ^ CXI at any prescribed value of f. 

This fact is related to an important difference between conformal geodesies 
and metric geodesies. It may happen that the parameter r on a conformal 
geodesic takes any value in R while the curve still acquires endpoints in M as 
T — > ±00. If the transformation in (R^ has a pole it may further happen that 
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the point sets spread out by x(r), x{f) coincide only partly, namely at the points 
where the transformation is defined, but each of the curves extends into a region 
not entered by the other one. In fact, there may occur an arbitrary number of 
such overlaps (we will encounter and make use of this phenomenon below). In 
certain contexts it may be preferable to call the union of the corresponding 
point sets a conformal geodesic. However, it will be more convenient for us to 
preserve this name for a solution of the conformal geodesic equation with its 
distinguished parameter. 

If the 1-form b is used to define a connection V along x{t) by requiring that 
its difference tensor with respect to V is given by 

V - V = S{b), (16) 

equation (|^) can be written in the form V^i — 0. Thus x(t) is an autoparallel 
with respect to V. Equation (^), which determines the connection V along that 
autoparallel, acquires some meaning by comparing it with (^. A congruence 
of conformal geodesies, covering smoothly and without caustics an open set U , 
defines a smooth 1-form field b and thus a connection V on U. If V, /, L 
are replaced in (^ by V, b, L respectively and transvect with x, we find that 
equation (^ takes the form of a restriction on the connection V in the direction 
of the congruence, given by L^i, x^ = 0. Here it turns out important that (^ 
does not contain the contraction of J^. 

Let / be an arbitrary 1-form and V = V + S{f ) the associated Weyl con- 
nection. Observing (^ with V, L replaced by V, L, we find for any curve x{t) 
in M and 1-form &(t) along x{t) the identities 

{VixY + S{b)x % = (Vii)^ -f S{b - /)a % x", (17) 
(Vi6),- i6^5(6)A^.i^-iA.i^ (18) 

= (Vi(6 - /)), -\{h- f)^ S{b - /)a % - ^A. i^ 

It follows that conformal geodesies are invariant under transitions to general 
Weyl connections and, in particular, under conformal rescalings of g in the 
following sense: if x{t), 6(t) is a solution of these equations with respect to V, 
then x{t), 6(t) — f\x{T) is a solution of the conformal geodesic equation with 
respect to V = V + S{f). Note that the parameter r is an invariant of the 
conformal class of g. It is determined by the initial conditions (|o|) but it does 
not depend on the Weyl connection and the metric in the conformal class chosen 
to write the conformal geodesies equations. 

Let Cfc, A: = 0, 1, ... n — 1, be a frame field which is parallel along x{t) for 
the connection V of ( p^ and satisfies g(ei, Cfc) = diag{l, ... 1, —1, . . . , —1) 
at a;(r*) with some O = 0* > 0. It follows then from (^), (||), that this relation 
is preserved along x{t) with the function Q = 8(t) satisfying 

Vie = e<fe,i;>, e(n) = e*. (19) 
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Consider again a congruence of conformal geodesies as above and initial data 
0* chosen such that the function Q satisfying (19) is smooth and positive on 
U. Then the associated 1-form is given in terms of the Levi-Civita connection 
V of the metric g — g, for which the frame is orthonormal, by 

/i(T) = 6(r)-e-ideU(,), (20) 

along x{t). Equations ( [l9| ) and ( pi] ) imply x(r) 

<h,x>=0, g{x,x) = Q'^ g{x,x) = g{x^,x^). (21) 

We assume the congruence, the function O, and the frame to be con- 
structed by a suitable choice of the initial data such that 

g{x,x)^l, eo = i. (22) 

From (13) and V = V + 3(0-^ dO) = V - S{h) follows V^efc = - < h,ek> 
X + g{x, Cfe) , where the index of h is raised with g. This implies that 

FxBk = ViCfe - g{x, ek) V^i + g{V±x, Ck) x 



= - < h,ek > X + g{x, e^) /i" - g{x, Cfe) + g{h^,ek) x = 0. 

Thus the frame ek, which is parallely propagated with respect toV, is in general 
not parallely but always Fermi-propagated with respect to V. 

The following general observation, which follows by a direct calculation, 
describes the sense in which Fermi-transport is conformally invariant. Let 9 > 
be some conformal factor and g, g metrics on M such that g = 9^ g. Denote by 
F, F the respective Fermi-transports. Let x{t) be any curve in M such that 
g{x,x) — 1. For any vector field X along the curve we then have Fgx{9X) = 
9^ F±X, where 9x is the tangent vector of the curve parametrized in terms of 
g-arc length. 

Applying this to the situation considered above, we see that the g-orthonormal 
frame Oe^ is Fermi-transported along the conformal geodesies, if the latter are 
parametrized in terms of g-arc length. 

Let f : M — > M be a conformal diffeomorphism of g such that f*g — VL^g 
with some function on M . Since conformal geodesies are invariants of the 
conformal structure, it is not surprising that / maps conformal geodesies into 
conformal geodesies. If if is a conformal Killing vector field its flow defines a 
1-parameter family of local conformal diffeomorphism and thus maps the set of 
conformal geodesies into itself. Though we will occasionally make use of this 
fact, it will not be demonstrated it here and we refer to |l^ for the formal 
argument. 



2.2 The conformal deviation equations 

Following the discussion of the Jacobi equation for metric geodesies, we derive 
now an analogous system of equations for conformal geodesies. Let x(r. A), 



8 



6(t, A) be a family of solutions to (||) , ^ depending smoothly on a parameter 
A. We denote the tangent vector field of the conformal geodesies, the deviation 
vector field of the congruence, and the deviation 1-form by 

X — drX = X, Z — d\x, B ~ \/z b, 

respectively. Considering x = x{t, A) as a map from some open subset of 
into M and denoting by Tx its tangent map, gives Tx{[dr, dx\) = 0, whence 
[X, Y]^VxZ~VzX = 0. This, (I), and (Vx - Vx) ^ ~ ^[x.z] Z = 
R {X, Z) X imply the conformal Jacobi equation 

Vx^xZ = R {X, Z)X- SiB; X,X)-2 S{b; X, Vx Z), (23) 

where {S{b;X,Y)Y = S{b) ^ XPV' for given 1-form b and vector fields X, 
Y . The last two terms on the right hand side of ( p^ ) indicate why conformal 
geodesies are potentially more useful than metric geodesies for the construction 
of coordinate systems. Under suitable circumstances the acceleration induced 
by the 1-form b may counteract curvature induced tendencies of the curves to 
develop caustics. 

Caustics of conformal geodesies can be more complicated than caustics of 
metric geodesies because for a given tangent vector there exists (essentially) a 
3-parameter family of conformal geodesies with the same tangent vector. More- 
over, to analyse them we need to complement the conformal Jacobi equation by 
a equation which governs the behaviour of B. Writing (Vx V z ~ Vx) b — 
'^[x.z] b = —bR{X, Z), we get from (^ the 1-form deviation equation 

VxB = -bR{X, Z) + {VzL){X, .)+L{VxZ, .) (24) 

+ i {B ■ S{b; X,.) + b- S{B; X,.) + b- 5(6; Vx Z, .)), 

where {B ■ S{b; X, .))^ = B^ S{b)p ^ ^ XP. We refer to the equations (H), (H) 
as to the system of conformal deviation equations. They form a linear system 
of ode's for Z and B along the curves x{t). 

2.3 Conformal geodesic equations on warped products 

For later applications we work out the simplifications of the conformal geodesic 
equations which are obtained in the case where g can be written as a warped 
product. Thus we assume that there exist coordinates x^ for which the set 
/ = {0, 1, . . . , n — 1} in which the index /i takes its values can be decomposed 
into two non-empty subsets Ii , I2 with I ~ Ii U I2, = /i ft /2 , such that in 
terms of the coordinates x"^, A £ Ii, and x"^, a £ I2, the metric takes the form 

g = g^u dx'^ dx" = Hab dx"^ dx^ + kcd dx" dx'^ , (25) 
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hAB = flABix^), kab = kab{x''), I < p = kabk"^ < U - I, f = f{x^) > 0. 

The Christoffel coefBcients are then given by 

f B ^ C = 2 ^^'^ {h,BD,C + hcD,B — flBCD), f (, ^ c = —h^'^ f f^D kbc, 
f B c = 0, fs " C" = 0, 
fB"c = fc°B = y /,B k"" c, f f, " c = - A:'"^ {kbd,c + kcd,b — kbc,d)- 

Those which cannot be derived by symmetry considerations from the coefBcients 
above, vanish. The curvature tensor has components 

R"" bcd\3] = R"" BCD[h], R" bcdIq] = c y Dd f, (26) 

R" bcdig] = R" bcd[k] - 2 /i^^ f,c f,D A;" [c kj^b, (27) 

where D denotes the /i-Levi-Civita connection with connection coefBcients Ta^ c 
Components which cannot be deduced by symmetry considerations from those 
above, vanish. In particular 

R^BCd\3]=o, R^bcd\g]=o, R''BCD[g]=o. (28) 

This imphes i^cf^] = 0) 

Lcd[g] = ^Rcd[k]- 
^ kcd R[k] - { ^ R[h] + ^ / DaD^J + DAfD^f } k,d, 

Lab\3] = I RAB[h] - ^ hAB R[h] - ^ { DaDb/ - ^ hAB DcD^f } 

-J^ { R[k] -P{P-1) Del D^f } hAB. 
The conformal geodesic equations take the form 

x^ + Tb^cx'^x^ + rb^cx^x^= (29) 

-2 {be x'^ + he x") + {hBC x"" x^ + fkbc x^) h^"" bo, 
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-2 {be + b, i=) + {hsc i"" + fhc i") ^ fc"" 6c, 

bA-Tc'' AX^bD-t,''Ai'bd= (31) 
(&c + be i") 6a - ^ {h''^ bebc + j^ k''^ bb be) Had + Lab [g] i'', 

ba-T,'' ai'bD-Tc\i^bd-rc\i'bd^ (32) 
(6c + 6, i^) 6a - i (/^^^ 6b 6c + 4 bb 6,) kac + iab[5] i'- 

3 Conformal geodesies on vacuum fields 

Though many of the subsequent discussions apply to more general situations, we 
shall assume from now on that dim M — A, sign[g) = (1, —1, —1, —1) and that 
g satisfies Einstein's vacuum field equations = 0. With the resulting sim- 
plification and the notation introduced above the conformal geodesic equations 



and the equation for the frame now take the form 

V±x + 2 <b,x> x-g{x,x)b^ =0, (33) 

Vi6- < 6,i > 6+ ig«(6,6)i^ = 0, (34) 

ViCfc-l- < 6, i > efc+ < 6; efe > i - g{x, Ck) 6" = 0. (35) 



Since these equations admit solutions with vanishing 1-form 6, metric geodesies 
form on vacuum space-times a subclass of conformal geodesies. 

3.1 Conformal geodesies on Minkowski space 

Denote by coordinates on Minkowski space (M, rj) in which the metric takes 
the form 77 = r]^^ dx^' dx" = [dx^f-idx^ f-idx^ f-idx^f and let =S^'k^^ 
be the orthonormal standard frame. There are various possibilities to deter- 
mine the solutions to ( |3^) , (|3^), (|35|). The framework of the normal conformal 
Cartan connection (cf.^H), in which conformal geodesies are defined as "hori- 
zontal" curves on a certain bundle allows a purely algebraic calculation of these 
curves, because the bundle admits in the case of the conformally compactified 
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Minkowski space a purely group theoretical description. Since the curves on the 
bundle space are known explicitly there only remains the task to calculate their 
projectioin onto the base space. We just give the result of the calculation. 

Besides the initial data ( p^ ) we prescribe the value 8* > and a Lorentz 
transformation which determines the initial value Ck* = Ai k Cj of the 
frame e^. The corresponding solution is given by 



e(r) = ( 1 + At < 6*,i* > +i (Ar)S(i*, i*) 6*) ) , (36) 



^''M=< + ^ (^Arii^ + i(Ar)S(i,,i,)6^:), (37) 



bu{T) = (1 + At < 6*, i* >) Ki, - ^ Ar ry'(6*,5*) i*^, (38) 



^'' = -^)^' A^)^J' (39) 

with a Lorentz transformation A{t) — (A^ k{T)) given by 

A{T)^A, + ATbixi 



-^(ATi, + -(Ar)2,^(i„i,)6«)(6,.A, + -Ar77«(6.,5.)it). 

Indices are moved in this subsection with the metric rj. 

Conformal geodesies which satisfy r]{x, i) = at one point coincide (as point 
sets) with null geodesies, those for which 6* = are Minkowski geodesies. If 
is space- or time-like we can assume, possibly after a reparametrization, that 
< >= 0. It follows that the remaining conformal geodesies fall into 

one of the following classes. If i* and bi generate a space-like 2-surface in the 
tangent space of a;*, the corresponding conformal geodesies is a metric circle in 
the plane tangent to that 2-surface. If and bt generate a time-like 2-surface 
and is space-like, the corresponding conformal geodesies is a space-like metric 
hyperbola in the plane tangent to that 2-surface. If is space-like and b^ is 
null, the conformal geodesies is a space-like curve in a null planed. 

We shall mainly be interested in the last case in which and bl generate a 
time-like 2-surface and is time-like. An example of such a conformal geodesic 
is given by the curve 

h\ < ^, 

\a\ 



<r) = (: 



+ 



a T 
2 



-,0,0), 
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which is a time-Uke hyperbola satisfying rj^i, x^{t) x^{t) — — a~^. 

We shall use now time-like conformal geodesies to construct coordinate sys- 
tems on Minkowski space where the parameter r on the curves defines a time 
coordinate while the other coordinates are obtained by dragging along spa- 
tial coordinates on S. Most important for us is the fact that the construction 
provides a conformal rescaling and extension of Minkowski space to which the 
coordinates are adapted in a natural way. 

We denote by the restriction of the standard radial coordinate r on 
Minkowski space to the hypersurfacc S = {x^ =0}. In spherical coordinates 9, 
(j) the inner metric induced on S then takes the form h = —{dri + rlda^) 
with da^ = dO^ -t- sin^ 9 dcj)^ the standard line element on the 2-sphere. If 
we write = tan(^) with < x < ""j find that the conformal factor 
6* = 2 (1 -f r1)^^ = 1-1- cosx allows us to realize the conformal embedding of 
S into the 3-sphere S* = 5 U {i} with line element 6^ /i = —{dx^ + siv? x^cr"^), 
where i denotes the point {x — '''} space-like infinity. 

To define a congruence of conformal geodesies and a conformal factor, we 
choose initial data on S as follows. At x* = (0, r^u) G S, u £ with |m| = 1, we 
set bl = au^ tc (sum over c = 1,2,3) and = 9;^^ eo such that ri{x, x) = 1. 
The function a on is determined by the following consideration. With the data 
above and r^, = we get 9(r) = 9* (1 — ^ 9~^) from (|^. To obtain, if 
possible, a 1-form b which is exact, we require 9~^ (i9 = & at a;* which gives 
a = 2r* {l + r^)~^. Our data are then spherically symmetric and we can express 
the conformal factor ( p6| ) and the curves ( |37| ) in terms of the time coordinate 
t — x'^ and the radial coordinate r to obtain 

To relate these expressions to known facts, we use r* — tan ^ and set ^ — tan | 
to obtain 

T 1 

Q = Quj, with n = cos s -I- cos X, uj = 1 + {- f ^ (41) 

2 cos^ I 

sins sinx 
, r= ^ . 42) 

cos X + cos s COS X + cos s 

Reading the last equation as a coordinate transformation, we get 

n^r]^n^ [dt^ - dr^ - da'^) = cje = ds^ - dx^ - sin^ x da^ . 

The set =M. x S endowed with the metric gE defines the Einstein cosmos. 
We see that ( p^ ) realizes the well known conformal embedding of Minkowski 
space into the Einstein cosmos, which maps the former onto the subset \s ± 
xl < TT, < X < of Me- The boundary dM of this set in Me supplies 
representations of future and past null infinity, future and past time-like infinity, 
and space-like infinity of Minkowski space, which are given by the subsets J^^ = 
{s±X = iTi"}, = {s = ±7r, X = 0}, i° = {.s = 0, x — ""I of dM respectively. 
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In terms of r and Q we get from (^) the metric 

Qh^ = ^^1^ dT^ - dx^ - sin2 xda^), (43) 

which extends smoothly to J^^ but does not extend to where uj — > oo. This 
is not due to an unfortunate choice of intial data on 5. The curve r — > z(t) = 
(s = 2 arctan^jX = 0) is a conformal geodesic on the Einstein cosmos which 
approaches as t — > ±oo. The freedom to perform parameter transformations, 



characterized by (14), (|15[), does not allow us to find a parametrization for which 
curve extends simultaneously to i~ and Thus the separation of i~ and 
realizes a conformal invariant. 

It is possible, however, to find reparametrizations under which the rescaled 
metric extends smoothly either to i+ or to i~ . The curve f z{f) = (s = 
s* + 2 arctan(^ — tan-^),x = 0) is again a conformal geodesic, because dg is 
a Killing vector field for g^- We choose s» to be a constant with < s» < 
TT. Then z{f) is related to z(t) by the parameter transformation t — f {1 + 
tan^ ^ — ■^T tan ^)~^- Comparing with ([l^) we are led to set = 0, = 0, 
Ae — (1 + tan^ and a* r|{x^,,x^,) = — tan ^ and to consider the initial 

data 

JL 1 - , ,^ ^ .1 1 - tan^ I 

^* = : — 2^^*' 0* = (?A2;*,a;*)2 = — 5-^6*, 6* = 5* - 



1 + tan^ ^ *' * *' * cos^ ^ *' * * 77(i*,i*) * 

(44) 



Observing these data and setting ^ = tan ^' + tan ^ , we get from {^7\ 
again equations (53), while (Bq) yields now 



tan % 1 „ tan^ 4f - tan^ ^ ^ 
^ l + tan^^ 4 (l + tan2^)2 ' 

with as in (|l]) and w = cos /"-^ . It follows that the metric 

e^rj^Cu^ ids' - dx' - rfa^) = Cu^ {\ df^ - dx^ - da^), 

extends regularly onto the domain |s — s*| < tt containing i+. 

If we insist on a conformal factor which defines a conformal compactifi- 
cation of {S,h) and initial data ensuring b = Q^^dQ on S, invariance under 
t —t implies the existence of a point on S where b vanishes. The conformal 
geodesic through such a point will then be a metric geodesic. In our first exam- 
ple such a point is given by = and it follows from ( |40| ) that the parameter 
T takes values in M and the conformal factor is constant on the curve r = 0. 
Thus, in order to achieve a compactification which extends smoothly to for a 
finite value of the parameter we have to choose non-time-symmetric initial data 
such as (^). 

Notice that could have been chosen to be a function on S. The possibility 
to change the parametrization while leaving the curves as point sets unchanged 
offers a large freedom to select slices of constant parameter value. 
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3.2 The conformal factor and the 1-form on general 
vacuum fields 



In the fohowing we shah derive some general, explicit information on the solu- 



tions to (|33|), (34), ( pSD (cf. the more general discussion in on solutions to 
the vacuum field equations i?^^ — Xg^^, with a cosmological constant A). 
From (H), (H) follows 

Vi < 6, i > = - < 6, i >2 +i g{x, x) g^{b, b), (45) 



V,~g\b,b)^<b,x> ~g\b,b), (46) 

where the index sharp on the symbol of a metric indicates here and in the 
following the contravariant version of that metric. 

Assume again that the congruence of conformal geodesies, 0, and the frame 
have been chosen such that (p2), whence g{x, x) — 8^^ holds along the congru- 



ence. Equation (11) then implies Vx Q ~ Q < b,x >. Taking a derivative and 
observing (psj ) yields V| 6 = 1/2 g*'(6, 6) 8^^. Taking another derivative and 
observing (Bm gives finally 

V| 8 = 0. (47) 

In terms of the initial data (|l^) this yields with 8, = y/ {g{x* , i,)) the explicit 
expression 

8(t) = 8* + At 8* + i (At)" 8* (48) 

= 8* ^1 + At < > -fi {ATfg{x*,x^) g''(6*,5*)^ . 

Similar steps lead to Vi {g{x,x) {g^{b,b))'^) = 0, whence 

Q-^g\b,b) = Q-^~g\b,,b,) = 2 8,. (49) 

In particular, the sign of g^{b,b) is preserved as long as 8 > 0. From (|^, 
(35) follows Vi(6 < 6, Cfc >) ~ l/2Qg'^{b,b)g{x,ek)- From this equation and 
(19), (^ ) the following explicit expression can be derived for the components 
bk = < b, ek > of the 1-form b in the frame 

ba = 8"^e, ba = O^^da with da = < 6*, 8, efe, > . (50) 

This implies with ( ^9| ) the relation 

8" - S"'' da db = 8" V'' bj bk = 8^ g^b, b) = g\b, 6) = 8 8:;^ g«(6,, 5,) = 2 8 8*, 

(51) 
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While the expressions (^8|) , depend on invariants built from the initial data, 
they are universal in the sense that their form does not depend on the solution 
9- 

Using the representation of the 1-form with respect to the connection of g, 
i.e. (pO|), we can write by ( pO[ ) on the congruence 6^(6) = 6<5, efc>— 0< 
/i, Cfc >. If there is a point on a given curve of the congruence where O vanishes 
and 6 and h remain bounded on the curve, (pO[), (51) imply 

V^ejCe) efc(e) ^0 as 6-^0. (52) 

3.3 The ^-adapted form of the conformal geodesic 
equation 

We write b = h + C,x^ (where here and below indices of 1-forms and vectors are 
moved with g) with h such that 

<6,i> = 0, whence C = g\h,b) = < b,x + g\b,b). (53) 

g{x,x) 



Equations (33), (|3J) are then equivalent to Vq± Qx — U, Vei 6 = —9^{b, b)Qx^ . 
We introduce the parameter transformation 

with inverse r = T(f) and set x(f) = x(T(f)). Then x' = dfX — Qx satisfies 
g{x',x') = 1 and we obtain the g-adapted conformal geodesic equations 

Vs'x'^P, V^,b^p^x'\ (55) 

where, by (|^), 

f3^ = -g^ (6,6)= S"'' da 4 = const. (56) 

along the conformal geodesic. These equations bring out the important role of 
6. If 6 vanishes at a point it vanishes along x{f) and the curve is a ^-geodesies. 
We determine the transformation (|5j) under the assumption that 

ff«(6,,6,) = .9(i,,i,)g«(6,,6,) <0. (57) 

It follows then from (|4^) that 0(t) vanishes at 

with T+ 7^ T_ , such that 

e = i .g«(6,, 6,)(t - t+) (t - r_), (59) 
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and the integration of ( |5^ ) gives 



nr)=n + ^log^^:i4(^, (60) 

1^1 [T- - T+) (r* - T-) 

whence 

. 2 9. sinh(f At) 

cosh(lfi Af) - < > sinh(l|i Af ) 

In terms of the parameter f we also get 

e = , ^ ^. (62) 

1 cosh(14i At) - e* < 6*, > sinh(14i Af )' 

Suppose that the solution admits a smooth conformal extension for which 
= on J''^, that the extension admits a point z+ such that J'^ coincides 
with the past light cone of i+, and that one of our conformal geodesies, x{t) 
say, passes through i+ for a finite value of r. Then condition ( p7[ ) precludes 
a dicussion of the zero of & on x{t). If ( |58| ) describes the zeros of 8 on the 
conformal geodesies covering a neighbourhood of x{t), then, approaching x{t), 
we should find t± Ti, and consequently (3 — 0, < b.^,,x^ > 7^ on x(r) and 
r- -T 2 

H — I* <b,,x,>- 

3.4 The ^-adapted conformal deviation equations 



Denote by x{t,X), &(t,A) a smooth family of solutions to ( |55| ) with family 
parameter A. As before, we write 

X = drx = x', Z = dxx, B = Vzb. (63) 

Following the derivation of the conformal deviation equations and observing that 
vacuum field equations = 0, we obtain the g-adapted conformal deviation 
equations 

\7xyxZ^C{X,Z)X + BK (64) 



WxB^-bC{X, Z) + {Wz9Hb,b))X^ +g\b,b)Wx Z\ (65) 

where C denotes the conformal Weyl tensor of g. 

While in general V2(g'*(6, 6)) 7^ 0, we have by ( |5^ ) always g'^(b,b) — const. 
along x(f), which implies Vx (VzgHb^b)) = Vz (Vx.9*'(6, &)) = 0. Thus the 
coefficients of X^ and Vx Z^ in the second line of ^ are constant and known 
along x{f) by their initial data at some f*. 
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3.5 Conformal geodesies and the (^-adapted conformal 

deviation equations on warped product vacuum fields 

It will be shown now that on warped product vacuum fields the discussion of the 
conformal geodesic equations reduces under suitable assumptions to the analysis 
of one equation and that a similar statement is true for the g-adapted conformal 
deviation equations. 

If g can be written in the form (p5|), the vacuum field equations take the 
form 

RAc[h] = j DaDcJ, Rac[k] = /"(P-') DAif^P'^'' D^f) kac, (66) 

and the dependence of the various fields in the second equation on a:° implies 

Rac[k] = ^ kac, R[k] = pf-^P-^^ DAif^"''^ f) = const. (67) 

We shall assume from now on that p = 2 and that the indices A, a take values 
0,1 and 2,3 respectively. Then RABCoih] = R[h]hA[c hD]B and RAc[h] = 
^R[h]hAc and thus 2DADBf = Hab DcD'-^ f by the first of equations ( pq) . 
Contracting with D^, commuting derivatives yields Dsif^ DaD'^J) = if ( |66[ ) 
is used again. It follows 



2c= p DaD^ f ^ const., DaDb! ^ -^Hab, ^ascdN = y| ^a[c 



(68) 



These equations imply with (p6|), (27) 



R BCD[g] = J^h [chD]B, R"' bcd[g] ^ — k" [ckd]b, R"^ BcD[g] = --p k" chsD- 

(69) 

It follows that i/Ap = unless c 7^ 0. 

If Einstein's vacuum field equation L^,^ = holds, equations (|29|), (^, 
(^, (|3^ ) admit solutions satisfying x°- = 0, 6c = 0, and these solutions obey 
equations which can be written in the form 

Dii=^-2<b,x> x + h{x,x)b\ Dib^ <b,±>b-^h^{b,b)x\ (70) 

where all quantities are derived from h. We shall discuss now the ^-adapted 
version of ([70|), assuming that A = det(/iAs) < 0. With 



eh^ y/\A\dx° Adx\ (71) 

it follows that e;i(i, .) = y/\A\{x^ dx^ — i^dx") and h^th{x, .),eh{x, •)) = 
—h{x,x). Since the vectors i, 5" are contained in the 2-dimensional space 
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spanned by Da, A = 0, 1, and < 6, i >= 0, the equations above and ( p6[ ) 
imply the representation 

h^±l3eh{Qx, .)=±l3eh{x', .), (72) 

where the sign is determined by the initial conditions and the convention for 
/3. Using this expression in the second of equations |5^, gives 0^ x' = Dx'b = 
± (3eh{Dx' x' , .), which can be rewritten in the form 

Ds'x' = bK (73) 

Thus, the g-adapted forms of the two equations ( [70| ) are equivalent to each other 
and we only have to solve ( |73| ) with b as in (|7^) to obtain a solution to ([70| ). 

Observing the formulae for the connection coefficients, the results (Um for 
the curvature tensor of warped products, and 

C^BCD[g]^^h^[chD]B, (74) 

equations (^J) , (|5|) are seen to be equivalent to each other and to the equation 

DxDxZ ^'^eUX,Z)eh{X, ± {Dz [3 ^hiX, + (3 eh{Dx Z, .)") • (75) 

The sign here is chosen as in ( fz^ ) and use has been made of the identities 
ehiehiX, .)"•) = XhiZ,X)-Z = eh{X,Z)eh{X, .)«. 

4 Conformal geodesies on the Schwarzschild 
space-time 

For our discussions various forms of the Schwarzschild line element will be 
needed. Its standard form with mass m is given for f > 2 m by 

il- —j dt^ - il- — \ df^-f^da^. (76) 

In terms of the retarded and advanced null coordinates 

w = i — (f + 2 m log(f — 2 m)), u — t + (f + 2 m log(f — 2 to)), (77) 

the line element (|7^) is obtained in the forms 

2 m \ / 2 TO \ 
1 ] dw'^ + 2dwdf- r'^da'^. 5=1 \ dv'^ - 2 dvdf - f'^da'^. 



r 



(78) 



respectively. These extend analytically into regions where r < 2 m. The re- 
tarded null coordinate w extends smoothly through J'+ and through the past 
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horizon to the past singularity, while the advanced null coordinate v extends 
smoothly through J^" and the future horizon to the future singularity. 
The transformation 

/ r \ ^/"^ f t / \ ^^"^ t 

s = I 1 exp(- — ) sinh(- — ), p — ( 1 expf- — ) coshf- — ' 

applied to (^6|) with m > produces the Schwarzschild-Kruskal line element 

32 m^ , f 



g = Kids'' -dp^)-f'^da^ with K = ^^^^ exp{ — ), r = k{p^-s'') 

r 2m 

(79) 



where k denotes the inverse of the map ]0,cxd[ 3 f — > — 1) exp(2^) G 
] — 1, oo[. This line element extends analytically to a non-degenerate line element 
on the domain p^ — > —1. 

The coordinate r, given for f > 2 m by 

r = — + — ^ resp. ^ ^ 2 ^ + •\/F(f--"2m)^ , (80) 
yields the isotropic Schwarzschild line element 



^ m 



2 



l + f) dt'~{l + ^ndr'+r'da'), r>^. (81) 

On the hypersurface {t = 0} it induces initial data which can be extended 
analytically to an initial data set 

(5 = {0}, h = -{l + ^ f {dr^ + da^), x = 0), (82) 

where x denotes the second fundamental form. These data may be identified 
isometrically with the initial data induced by the Schwarzschild-Kruskal metric 
[79| on the hypersurface {s — 0} by the transformation 

.1 )■ 



\/2mr 4mr 

In the following discussions one may think of the conformal geodesies as be- 
ing constructed on the Schwarzschild-Kerr solution ( [79|) with initial data being 
prescribed on the hypersurface S* = {s = 0}. It will be convenient, however, 
to specify them in terms of the coordinate r in (|8|). To discuss the different 
aspects of these curves we will use those coordinates which will appear to give 
the simplest formulae. Because of the symmetries s — > —s and p —p, it 
will be sufficient to consider the region {s > 0, p > 0} if we prescibe data for 
the conformal geodesies which respect these symmetries. The data will also be 
spherically symmetric. 

The function f extends by the first of equations ( pO[ ) to an analytic function 
on S which takes its minimum fmh-i = 2m at the "throat" {r = rh = The 
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space {S, h) has two asymptotically flat ends. In terms of ( |8^ ) the end it at 
r = oo has the "usual" representation, while the end 12 at r = has a coordinate 
representation adapted to the metric {1 -\- ^)~^ h which may be thought of as 
a conformal compactification of h at the end 12- The map 

,TO,9 1 , , 

''-(y)-' (83) 

which corresponds to the isometry p —p allows us to show the equivalence 
of the spatial infinities ii and «2- It leaves r invariant and has {r = r/j} as its 
fixed point set. 

4.1 Conformally compactified Schwarzschild-Kruskal data 
and initial data for a congruence of conformal geodesies 

We shall in the following prescribe initial data on S for a congruence of conformal 
geodesies whose ^-adapted initial vectors x' coincide with the future directed 
unit normals of S. By f* will be denoted the restriction to S of the function f on 
the Schwarzschild-Kruskal space-time and r will be considered as a coordinate 
on S. For the conformal factor we choose 

which implies 

-elh^ {r+'^y^idr^ +r^da^). (85) 

It follows that the transformation r — tan ^ , < x < tt (and, to get simple 
equations, ^ — tan^) could be used to make the conformal compactification 
achieved by 8* manifest in terms of coordinates, since it realizes an embedding 
of iS* into the unit 3-sphere with the poles at x = 0, x = tt corresponding to 
the ends 12, ii respectively. 

The choice (|4|) is particularly well adapted to the Schwarzschild-Kruskal 
geometry and implies simple formulae. However, it is important to note that 
our basic results do not depend on it. For the analysis of the field near space- like 
infinity other choices might be preferable (cf. [^), because the metric (|8|) does 
not extend smoothly to space-like infinity. 

We set furthermore 

2 fr — — ) 

h = h = e:Ue, = - ^l{ dr, (86) 

and, observing F = 1 — 2?7i/f* = (r — + choose 

2r(r-f) 
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such that (3 is analytic on S and [3 = \J —g^(b^,, 6*) = -p- ■\/F(r*) > near the 

end «i in analogy to our procedure on Minkowski space. If h is used to raise 
indices, b\ is outward pointing at the ends ii, 12, as it should. With the choices 
above the general formula (|59| ) reduces to 

Q^F{f.) (^(^f-r^), (88) 

where functions with subscript * are assumed to be constant along the conformal 
geodesies. Note that our initial data and gauge conditions are preserved by 
the isometry (^3|), and are adapted to the spherical symmetry and the time 
reflection symmetry. We have (3 — > —(3 under (|83| ) and /5(r) = precisely for 
r = Y • This means that conformal geodesic satisfying the initial data above at 
points with r = ^ '^iH be metric geodesies in the hypersurface {p — 0} of the 
Schwarzschild-Kruskal space-time. 



4.2 The conformal geodesic equations on the Schwarzschild 
space-time 

To discuss conformal geodesies on the Schwarzschild space-time we write the 
g-adapted form of the conformal geodesic equation in terms of the line element 
([76|), which can be written in the form (p5| ) with the metric h given by 

h^Fdt^~—df^, F = (I —). 

F r 

Initial data for the conformal geodesies will be given on the hypersurface S = 
{t = 0} with r, = 0, f, = on S. 

Because of ( [72|) (where due to our conventions we have to use the minus 
sign), equations (|73|) take the form 

t" + E^r't'^^Pf\ (89) 



f"-^{f'f + ^{t'?-Fpt', (90) 



with (3 satisfying (p6D. Assuming on S that the initial vector is the future 
directed unit normal to 5, the initial data on S are given by 

t* — 0, f* > 2m, = _ , = 0, 6t, — 0, 6f* — —f3{f^) 



The ^-normalization gives 

Fit'r-^{fr = i- (91) 
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Solving for i' > and inserting the result into equation (pG) gives 



= f" + - F,r -p y/F+if'y. (92) 



Dividing by the square root and multiplying with r' leads to d{y/ F + [r'Y 
(3r)/df = 0, whence 



y^F+{f')^-pf = -f, (93) 



with the constant of integration given by 7 = ^/Fir^ — P{r^) r*. It follows that 



f' = ± V(7 + /3f)2 -F(f), (94) 

with a sign which depends on the value of . 

Given f(f), we obtain t{f) by integrating {t')^ — (fif + ^)'^ F^"^ for the given 
initial conditions. However, the integration of t is in general not particularly 
interesting, because t does neither extend smoothly through nor through 
the future horizon. 

To study the extension of the conformal geodesies through it is conve- 
nient to use the first of the line elements ( ^8| ) . The conformal geodesic equations 
then read 

w" - ]^F,f {w'f = -f3w\ r" ^i^FF,f {w'f+F,f f'w' = (3{f' + Fw'). 

These equations imply for r the equations obtained above. The normalization 
of the tangent vector reads F (w')^ + 2w' f' — 1. Solving for w' and requiring 
w' > on S gives 

which has to be integrated with the initial conditions = — —2m log(f* — 
2 m). Thus w{f) can be obtained by a simple integration once f(f) has been 
determined. 

To study the extension through the horizon, it is convenient to use the second 
of the line elements (^8|) . The conformal geodesic equations then read 

v" + i F,r {v'f =(3v', f" + \FF,r {v'f - F,r f'v' = -(3 (f' -Fv'). 

The normalization of the tangent vector gives now F {v')'^ — 2v' f' ^ 1, which 
leads to the equation 

1 



v' = j{VF+{f'r + r'), (96) 
which has to be integrated with the initial conditions = f* + 2 m log(f* —2 m). 
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4.2.1 Conformal geodesies on which r is constant 



The expUcit solution of the key equation (94) requires the discussion of three 
different cases. We first consider the borderhne solution. 

A change of the sign in (H) should occur near conformal geodesies along 
which f is constant. Requiring f' = in (^2|) gives 1/2 F' = P^/F. With /3 



given by (87) this condition has the solution 



3± ^/5 



r = — m resp. r ^ r+ = m. 

2 4 



Inserting r' = into (|9l|), using (^) with = , = , and observing (| 
( p7| ) with r = r-|_ gives on the conformal geodesies on which f = f 

f 25 , /2e*+/3T 

t — — ^ = — m log 



and thus i — > oo as r — > r,- = ^ ^7^^!'* . 

It follows now from equations ( |92| ) , (|9J) that each of the conformal geodesies 
specified by our data falls into one of the following four classes: (i) the confor- 
mal geodesies passing through points with r = which coincide with metric 
geodesies in the Schwarzschild-Kruskal space-time and approach the singular- 
ity, (ii) the conformal geodesies passing through points with r = r±^ which are 
tangent to the static Killing field dt and approach time-like infinity for the finite 
value Ti of their parameters, (iii) the conformal geodesies passing through points 
with < r < r_ or with r+ < r < cx), for which r is monotonically increasing, 
(iv) the conformal geodesies passing through points with < r < ^ or with 
Y < r < r+, for which f is monotonically decreasing. Though it will be seen 
below that the integration procedure for ( p4|) depends on the classes specified 
above, it is clear that our data, which are smooth on S, determine a smooth 
congruence of conformal geodesies which is free of caustics near S. 

4.2.2 Conformal geodesies on which f is increasing 

Because of the symmetry of the data it is sufficient to discuss the case r+ < 
r < oo. Under this assumption the "+" sign holds in (|9j). The radicand in this 
equation factorizes as 

ij + f3ff ~F{f) = y (f-r*)(f-a+)(f-a„). (97) 



Because 7 = — ■\/F(r*), it follows that 



a± — ±a with a = 



2F{f^y 



Since a < f.^, and a — > f as ^ r, the polynomial on the right hand side of 
( p7[ ) has under our assumption three different zeros while in the limit above two 
zeros will coincide and the integration procedure needs to be changed. 
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The use of the notation 

{a + n)f) ' ~ \2 
in (|9j) gives, after an integration, 



pf=- V(e2 - 1) (e2 - fc2) 



i+f 



(1 - 62^2) ^(l-t2)(l_fc2 



where, by our assumptions, 1 < e < ^/2, < k < 1. It is well known how to 
express the elliptic integral of the third kind which occurs above in terms of theta 
functions (for all manipulations involving elliptic integrals, elliptic functions, 
theta functions, etc. we refer to ^). The substition t = sn{u, k) with Jacobi's 
elliptic function sn yields 

u 

dv 



/3f^-V(e2-l)(e2-fc2) / - (98) 



K 



e y (1 — e2 sv?v) 



i(2T?("-c)) 



= -2 Z{c) {u - K) + log 
Here Z denotes Jacobi's zeta function 

Z[u) = - (log 9,{- -)) = - sm(. -) )^ i_2,2„-i eos(.f) + g4»-2 - 

Thus Z is an analytic function on the real line with Z{u) > for < u < 
K and Z{0) = Z{K) = 0. The theta functions are given in their product 
representations by 

oo 

ei{z) = 2qi sin z J| (1 - ») (i _ 2 g2« cos 2z + q^"), 

n=l 

oo 

02{z) = 2qi cosz [](l-g2")(l + 2g2»cos2z + g4"). 



04(Z)= J|(l-g2„^^^_2^2„-l cos2z + g4"-2)^ 
n=l 

with q — exp(— TT ■^). The constant K > is the complete integral of the first 
kind 



K{k) 



dt 



^(1 - t2) (1 - p ^2)' 
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K' = K{k') with k"^ = 1 — fc^, and c is the (unique) number in the open interval 
]0, K[ satisfying sn (c, k) = 1/e. 

With n = 0, = 0, and @ we get from equation (|o|) 



/3f = -log.9(T) with g{T)= 'i^*^J . (99) 

2 H* + p T 



Setting now a; = i^T — m to exhibit the symmetries of the solution, we obtain 
T = G{x,k) 



2 9, 02(^^)e"^(^)-02(^^)e-"^(^) 



/3 02(^i^)e^2(c) +02(i^)e-xZ(c)■ 
Since the denominator of the function G is positive in an open neighbourhood 
of the closed interval [~{K — c), X — c], the function G is analytic there. While 
T ^ oo as M — > c and r — s- — oo as u ^ 2 iiT — c by (98), we have now r ±^-^, 
whence f ^ ±oo, as a; ^ ±{K — c). It follows from the ODE satisfied by f 
that we can solve the equation above for x(r) with re] — ^^[. A direct 
calculation shows that G'(a;, A:) > at cc = ±(/f — c). This implies that x{t) 
extends as an analytic function into an open neighbourhood of 

Inserting now x{t) into the equation t = sn{K ~ x,k), we finally get 

n_ 2(1-P)fc^ sn'ixjr)) 

f 2fc2-l l-fc2sn2(x(r))' ^ 

In the limit when r ^ r+ we have fc — > 1 and thus sn (u, k ) — > tanhu. The 
formula (IOC) thus suggests that the solution f{T,r) of (|lO0| ) approaches the 
constant solution r = f in that limit, as we know already by general arguments. 
However, since K ^ oo while c 1/2 log((\/2 + 1)7(^2 - 1)) as A: -> 1, the 
precise behaviour of the solution in that limit requires a quite careful discussion 
involving also the behavio ur of G(x, k). 

The right hand side of ( |lO0| ) is an analytic function of r in an open intervall 



containing [~Tscriir),Tscri{r)], with Tscri{r) = ^ = (^r+^){r-^) - vanishes 
precisely at the points r — ±Tscri at which the conformal factor ( ^8|) vanishes. 
Furthermore, it depends analytically on r for < r < oo. It follows that 

f(r, t) 8(r, t) is positive and analytic for r G]r_|_,cx3[, t G [—Tscri{r),Tscri{r)]- 

(101) 



Using the coordinate z = ^ in the first of the line elements (78) and rescaling 
with the conformal factor = z gives the smooth conformal representation 

g = {l~2mz)dw^ - 2dwdz - da^, (102) 

of the Schwarzschild metric which extends analytically through future null in- 
finity, given here by J^'^ ~ {z ~ 0}. 

Integrating (|9^) with f given by ( |10C| ) and writing the solution in terms of 
z, we obtain our conformal geodesies in the form 

r— >(w(r, r), z(t, r)). (103) 
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From the discussion above it is clear that the first of the functions on the 
right hand side is an analytic function of both variables for r €]r+,oo[, t S 
[—Tscri{r) , Tscri{r)]. Wc show that this holds true also for the second function. 
It is clear that the function w is analytic near S. Parametrizing w in terms of 
z, we obtain from (p^), ( p5| ) 

4^ - - (1 - n z) (1 - z^) + (/? + 7 2) (1 - z) (1 - a2 ^2)^ . 

The assertion now follows because the function on the right hand side is analytic 
for z in an open interval of the form ] — e, with some e > 0. 

It follows in particular that there exists a smooth function w{r), with r g 
]r+,oo[, such that w{T,r) w(r) as t ^ Tgcri (or, for symmetry reasons, as 
T —Tscri)- We show that 

w{r) ^ oo as r ^ r^,, w{r) — > —oo as r — > oo. 

The first assertion follows from the observation that 

w = —r — Ira log(r — 2 m) H m log — — — 

4 \2 0* — p T 

along the conformal geodesic with r — f and the fact that the solutions are 
jointly smooth in the initial data and the parameter. The second assertion 
follows from a comparison of w with the solution to 



which satisfies u — w aX z — Since > 4^ > 4^ for z e [0, ^1, it follows 
that the value u{r) of m at z = gives an upper estimate for w{r). Since the 
direct integration gives u{r) with u(r) — > — oo as r ^ oo, the assertion follows. 

4.2.3 Conformal geodesies on which f is decreasing 

It is sufficient to discuss the case ^ < < Now ( p^ ) must hold with the 
"— " sign. The function a in the factorization ( p7| ) then* satisfies f * < a and 
a — + oo as r* ~> Y while a — > f as — > f . If we set 

such that e>l,0<A;<l, the integration of ( p^ ) yields 
/3f = v/2(2-fc2) 



i+f 



(^2 e2 - 1) ^(t2 fc2 _ 1) (^2 „ 1) 
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The substitution k sn{u, k) — 1/t and c €]0, K[ such that sn{c, k) — 1/e give 



K 

2 — k^ f sn?v 



^^dv^-2Z{c){u-K)+\og 



2 



2^4" 



Setting X = K — u and observing (pa), we finally obtain 



T = F x, A: = ^ tanh k^ \ — — / —- —dv . (104 

^ ' P \ 72-fc2(l + sn^v) ' ^ ' 

The function F{x) is real analytic on M with F'(a:) > and F' [x) = iff 
X = (2 m + 1) if with to G Z. Thus (104) can be solved to obtain an analytic 
function x(t) which maps the interval ] — Ts,rs[, with 2Qs,/(3 < Tg = G{K,k), 
diffeomorphically onto ] — K,K[. The solution can then be written 

, _ , (2-P)(l-.n^(x(r))) 

* 2 ~ k^l + sn^xiT))) ■ ^ ^ 

From the discusssion above it follows that f(T) and df/dr =Foo as 
T — » ±Ts. Thus in this case there does not exists a smooth (though a continuous) 
extension of f(r) beyond its physical domain. The limit r — > r+ implies fc — s- 1. 
In particular, it follows from ( |104| ), (105) that the right hand side of ( |l05| ) goes 
in this limit to f for constant r and Tg ^ Ti. 

The limit ^ implies k ~* Q. Since (3 — Q aX r — the expression 
for Tg appears to give a nonsensical result at that point. However, we have 
2G* fc^//3 2/m. Expanding the right hand side of ( |104| ) and observing that 
K{0) = tt/2 thus gives Ts it /Am as r ^ "i/2, consistent with the fact the 
conformal geodesies approach in the limit metric geodesies with length f = 

„ „ — 7 = TO TT. 

0.^(r*— 2m) 

To follow the conformal geodesies through the horizon, one has to integrate 
for given fir) the equation 



, l + l3r-^{-1 + l3ff~F{r) 

V = JTTTT , (106) 

F{r) 

whose right hand side defines a positive analytic function of f for f > 0. The 
conformal geodesies are then obtained in the form 

T ^ x(t, r) = (z;(r, r), f(T, r)), (107) 

where the function on the right hand side are analytic for (r, r) with r e] ^ , ^+ [, 
re [0,n.(r)[. 



4.3 Analytic coordinates covering the Schwarzschild- 
Kruskal space-time and its null infinity 

If we set now — r, y'^ = (j>, ~ 9 on S, drag these coordinates along with 
the congruence of conformal geodesies constructed in the previous section, and 
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set = T, we obtain a smooth coordinate system near S. The purpose of the 
following discussion is to establish the following result (we ignore here the co- 
ordinate singularity arising from the use of coordinates on 5*^ which can easily 
be removed): 

The conformal Gauss system defines a smooth global coordinate system 
on the Schwarzschild-Kruskal space-time which extends smoothly to null infinity. 



It remains to be shown that the conformal geodesies of the congruence do 
not develop caustics. Because of the symmetry of the functions and of our 
family of curves under s — > —s, p — *■ ~p, it is sufficient, to control the behaviour 
of the congruence on the subset {s > 0,p > 0} of the Schwarzschild-Kruskal 
space-time. Since the coordinates are regular near 5*, they are known to be 
regular in particular near the set {s = 0,/5 = 0}. Therefore we can work with 
the coordinates f or z, f or v, f respectively. 



Using the functions on the right hand sides of (103) and (107), we can in 
principle express the corresponding line elements in terms of the coordinates 



Though the resulting metric coefficients will be smooth where (103) and ( |107| ) 
are analytic, the Jacobian of the transformation may drop rank at various places 
and we may end up with a degenerate representation of the metric. Having 
available the explicit solution, we could try to check this by a direct calculation. 
Since the explicit calculation will be tedious and, in particular, because this 
method will not apply to more general cases, we prefer to employ an argument 
which is similar to the analysis of the behaviour of metric geodesies congruences 
in terms of the Jacobi equation. 

Consider the transformation provided by ( |107D . We need to show that the 
tangent vector field x — x! = O^^ X and the connecting vector field x^r = 
= Z — xT^r are linearly independent on their domain of definition. In 
terms of the 2-form ( [7l| ) (with h = Fdv^ — 2dvdf,x^ = v,x^ — f) this can be 
expressed as the requirement that the invariant e^iX, Z) does not vanish. 



In the case of (103) the nondegeneracy up to and in fact beyond null infinity 



would follow if it could be shown that 

~{fQ)'^{wZ^r-ZW^r)^0 for e]r+,Oo[, T € [0,Tscri{r)]- 

If z is replaced by f again, this condition translates in terms of x(t, r) = 

(u(t, r), f(r, r)) and ( |7l| ) (with h = F dw^ + 2dwdr, — w, x^ = f) into the 
requirement that 

^ £6^/1(2;, X^r) = 6 ^hiX, Z), 

in the domain given above. Here the factor Q is of course most significant 
because it vanishes at Tgcri- 

The proof that these requirements are met by our transformations relies on a 
differential equation satisfied by eh{X, Z). To make use of (|75|), we observe that 
the various representations of the Schwarzschild metric with warping function 
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f = r and second metric k = —dcr^ give the value c = i DaD^J = — m for 



2 

the constant in (1391). It follows 



Dx eh{X, Z) = eh{Dx X, Z) + eh{X, Dx Z) = e^(-/3e(X, .)", Z) + eh{X, Dx Z) 
= -l3h{X, Z) + eh{X, Dx Z), 

and similarly 

Dx Dx eh{X,Z) = 
eniX, Z)-2f3 h{X, Dx Z) + Dx Dx Z) 

= {P"" + ^)eh{X,Z) + Dzl3, 

where in the second equation Dx Z = Dz X, h{X, X) = 1, and ( |75| ) have been 
used. The invariant efi{X,Z) thus satisfies the ODE 

Dx Dx euiX, Z) - {(3^ + k) en{X, Z) ^ Dz /?, 

with k = ^IP^, and on S the initial conditions 

eh{X,Z), = ^ 5^>1, {Dxth{X,Z)),^Q for r > -. 

The quantity DzP, which is constant along the conformal geodesies, is given by 

DzP = I3,r = -2 ■ 



(r+f)4 • 

It vanishes at r± = tti, where f{r±) — 3m > and satisfies 

771 

Dz/? > for Y<r<r+, Dz/3 < for r+ < r. 

A lower estimate for eh{X, Z) can be obtained as follows. Denote by u and 
V the solutions to the ODE problems 

w" - iP^ + k)w ^ f, w(0) = 1, w'(0) = 0, 

where / = in the case of u and / = — 1 in the case of v. Then u is strictly 
increasing with u > cosh(/5f) for f > 0, and eh{X, Z) can be given in the form 

eh{X, Z) = u (euiX, Z\ + (1 - ^) Ax) • (108) 

Since [u — v)" — {0^ + k) {u~ v) = 1 and the fmiction u — v has vanishing initial 
conditions at f = 0, it follows that u > v for f > 0. Since v changes its sign, a 
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further estimate is needed. We derive a representation of v in terms of u. Since 
M > there exists a function / with v — f u. The ODE's satisfied by u and v 
imply for / the equation 

/" = -2-/'--, 
u u 

which has because of the initial conditions for u and v the solution 

Jo w Jo 
Since u is strictly increasing, it follows 

2/ T'e-'^^'dr' = -^(1- (/3f + l)e-^^) for f > 0, 
Jo P 



< 

which implies 

V 2 
- < — 

u- f3' 

A direct calculation gives —1 < 2/3^^ Dz /? < for r > r+. It follows that 



0<l--<^. (109) 



V Til 

e^{X, Z), + {!--) 13^^ >e,,{X,Z), for -<r<r+, 



e,,(X,Z), + (l--)/3,^>e„(X,Z), -1 = for r+ < 



Since ( |62[ ) gives under our assumptions 

cosl7(f f) 



,2%-. ^ (110) 



(108) implies that for given r > ^ there is a constant c > such that 

ee.(x,z)>e.cii^^^>e.c. 

cosh (^ T) 

In the region covered by ( [107] ) it suffices of course to get a lower estimate for 
eh{X^Z), because Q is positive where the conformal geodesies approach the 
singu larit y. On the curves with r = f, which seperate the domains where (103) 
and (107) are valid, (? = 0^ + k with c — const. > and thus 

eehiX,Z) = e{eh{X,Z)^ cosh{cf) + c-^ Dz [3 {cosh{cf) - 1) 

> e,ehiX,Z), > 0. 

Since P — but k > k^, = (2to)^^, Dz P = on the curves starting with 
r = ^, it follows that O = G* = (2m)^^ and a similar result is obtained. 
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4.4 Global numerical evolution for a class of standard 
Cauchy data 

The methods described above offer a possibihty to study in detail the complete 
numerical evolution of three-dimensional space-times without symmetries which 
are determined by certain asymptotically flat Cauchy data. In has been 
shown the existence of smooth, time-symmetric, asymptotically flat solutions of 
the vacuum constraint which coincide with certain given time-symmetric initial 
data on compact sets and with Schwarzschild data in a neighbourhood of space- 
like infinity. If these data can be constructed numerically, it is easy to determine 
numerically hyperboloidal initial data implied by the Cauchy data. 

Consider the time symmetric initial data set (p2[). The set S can be em- 



bedded by the transformation r = tan ^ into the 3-dimensional standard unit 
sphere S = S^. With the conformal factor Gf = 2 (1 r^)"! (1 + fi)-^ = 
2 cos^ f (1 + ■fcot^)^'^, the rescaled metric (9f )^ = -{dx^ + sin^x^f^) 
coincides with the restriction of the standard metric on to the set 5\ {io, ^ir} 
, where io = {x = 0}, = {x = tt}. If we set m = 0, we get Minkowski data 
and conformally compactified Minkowski data respectively. 

For given xo with ^ < xo < let ^ be a smooth function on R such that 
£,{x) = for a; < ^, ^(x) = 1 for x > xo and such that ^' > 0. We define V' G 
C°°(] — oo, 7r[) by setting it equal to 1 for x < -I and equal to 1 -t-^(a;) (cot § — 1) 
for I < a; < TT. Then V'' < and i < M = sup^<^ < oo. 

Suppose we are given time symmetric initial data on M.^ which agree with 
Schwarzschild data of mass m < jj near space-like infinity and are such that 
the metric h, suitably written in terms of the coordinates 0, 9, x on S, satisfies 

= (dx' + sin^Xtia^) for x > Xo- (HI) 

4 cos* ^ 

The conformal factor 

2 cos2 I ^ ^ 

is smooth on 5\{i7r}, coincides for x ^ Xo with 0f, has non- vanishing gradient 
in 5* \ {iQ, ijr} and goes to zero at i^. It defines a conformal compactification of 
the data such that h = Q^h coincides with the standard metric of the 3-sphere 
for X > Xo- We choose initial data 6* for the 1-form field which annihilate the 
normals of the initial hypersurface and satisfy 6* = dO* in 5* \ {tTi}- 

Since the time evolution of the data will be Schwarzschild near it can 
be determined there by the methods described above. It is clear that we can 
construct a smooth hyperboloidal hypersurface, which coincides with S on the 
the set {x < Xo} and extends to the future null infinity of the Schwarzschild 
part. It should not be difficult to determine the corresponding initial data for 
the conformal field equations, possibly by a numerical integration (as shown in 
[||, this reduces to solving a system of ODE's). Since there are codes available 
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to evolve such data numerically (cf. |8| and the references given there) we 
could in principle calculate their evolution in time. 

If such data, satisfying for a fixed xo, can be constructed for sufficiently 
small mass m they will be close to corresponding Minkowskian data and the 
result of suggests that they will evolve into solutions possessing complete 
past and future null infinities and regular points at future and past time-like 
infinity. If we use a gauge in which picks up on 5'\ {ijr} a suitable component 
in the direction of the tangent vectors of the conformal geodesies, it will be 
possible to construct a gauge of the type considered above which smoothly 
covers the complete future of the initial hypersurface as well as J""*" U The 
work in ||^ suggests that such solutions can be calculated numerically in their 
entirety. 

The numerical calculation of such space-time is, of course, not an end in 
itself. In fact, the solution will have somewhat curious features. It follows from 
[H that they will have vanishing Newman-Penrose constants. If the solution 
admit regular points it then follows from that the rescaled conformal Weyl 
tensor will vanish at those points. This situation, which is more special than the 
one considered in , suggests that the method of gluing a Schwarzschild end to 
a given solution of the constraints produces data of rather restricted radiation 
content. However, such calculations will allow us to perform detailed tests of 
the code under completely controlled assumptions and to study the robustness 
of the code and of the gauge conditions by choosing data with an increasing 
value of the mass which eventually yields the time evolution of a collapsing 
gravitational field. 

5 Concluding remarks 

We have described in detail the construction of a global system of conformal 
Gauss coordinates on the Schwarzschild-Kruskal solution which extend smoothly 
and without degeneracy through null infinity. Furthermore, we have shown that 
the conformal factor naturally associated with this system defines a smooth 
conformal extension of the Schwarzschild-Kruskal space-time which gives to null 
infinity a finite location in the new coordinates which is determined by our choice 
of initial data. 

We did not try to work out in detail the behaviour of the fields in the 
conformal extension constructed here. An analysis of the fields near space-like 
infinity can be found in H]. The behaviour of the fields near time- like infinity, 
which is of considerable interest for the numerical calculation of space-times, 
has still to be investigated. 

There is a property of the Gauss system which we only indicated but which 
may turn out to be quite important. While the regularity of a conformal Gauss 
systems is essentially decided by its underlying conformal geodesies (considered 
here as point sets), there always exists a huge class of different time slicings 
based on one and the same underlying congruence of conformal geodesies. The 
consequences of this freedom still have to be explored. 
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Of course, there are many more such coordinate systems. It would be in- 
teresting to see whether the initial data for congruences of conformal geodesies 
which lead to such coordinates can be characterized in a general way. 
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